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Abstract 

We show how to construct 5*0(4) invariant functions for strongly correlated Fermi 
systems on lattices of finite sizes. We illustrate the method on the case of the ID 
Hubbard chain with four and with six sites. 
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During the past years a great deal of interest has been devoted to the study of strongly 
correlated Fermi systems because of their possible role as models of high T c superconductivity [Q. 
Among these systems, the Hubbard model is certainly the simplest non trivial model for interacting 
electrons in a solid. In spite of its apparent simplicity, the mathematical and physical properties 
of this model, in dimensions higher than one, are still poorly understood. On the other hand, the 
huge dimension of the Hilbert space (4/ for a lattice of / sites) severely restricts numerical studies 
to clusters of small sizes. In order to reduce the dimension of this space one may use the symmetry 
properties of the Hamiltonian. Basis functions which account for the conservation of the number 
operator, for the conservation of the z component of the total spin S and for the translational 
invariance, were constructed by many authors ||. For the Hubbard model on bipartite lattices, 
however, one would like to use as basis functions the simultaneous eigenfunctions of the SO (4) 
(spin S and pseudospin J) algebra as well as of the translation operator. 

The aim of the present paper is to show how to construct SO (4) invariant functions for strongly 
correlated Fermi systems on arbitrary lattices. As a result we give a set of rules for constructing 
such functions which are easy to implement on a computer in terms of symbolic languages. We 
illustrate the method on the example of the Hubbard chain with / = 4 and / = 6 sites. For 
/ = 4 we find that the SO(4) diagonalization decomposes the original 256 * 256 hamiltonian 
matrix into blocks of dimension less or equal to four. The ground state at half filling is obtained 
by diagonalizing a 3 x 3 block corresponding to S = J = and with momentum equal to k = tt. In 
the case of the / = 6 chain we find that the ground state at half filling is obtained by diagonalizing 
a 14 * 14 band matrix corresponding to the S = 0, J = 0, k = tt space, with the width of the 
band equal to six. Finally, at the end of the paper, we estimate the sizes of the SO (4) blocks as a 
function of /, J, S. 

Let us start by introducing the Hubbard Hamiltonian as 

/ / 

cri i 

where cj , Cj a , (a =f or J, ), are usual fermion creation and annihilation operators. As well known, 
this model is invariant under SU(2) spin rotations with generators 

f i 
5+ = E c n c n' S- = (S+)K 5 z = -(iV T -iV;). (2) 



and under SU (2) pseudo-spin rotations with generator 



j + = T,(- 1 H c t' j- = E(- 1 ) i ^^ J ^^- (3) 

i i 

One readily checks that S, S z , J, J z , H form a complete set of compatible observables. Moreover, 
from the expression of S z , J z in Eq.s( |^,^) it is evident that the states of the system must have 
spin and pseudospin both integer or both half integer, this giving a Z 2 reduction which leads 
to a global 5*0(4) symmetry [Q. To construct the basis functions which span the irreducible 
representations of SO(4) we will use the representation theory of the permutation group. As well 
known, Young tableaux are natural objects for this purpose. To this end we remark that the 
commutation between the permutation operations and the spin rotations implies that the irreps 
of Sf are automatically irreps of the SU(2) spin algebra. This was used in ref. |^|, ||) to solve the 
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Hubbard model with unconstrained hopping for clusters of arbitrary sizes. Here we extend this 
result to the SO (A) symmetry and to arbitrary lattices. To this end let us denote with the symbols 
|3 >, |2 >, |1 >, |0 > the four states on a given site (respectively the doubly occupied, the single 
occupied spin up and spin down states and the vacuum) and let us introduce the quantum number 

M=~(f-2J Z ) + S z . (4) 

Eigenmanifolds of S 2 ,S Z ,J 2 , J z can be constructed from the irreps of Sf by considering all possible 
partitions (mi, TO2, m/) of M into / parts (compatible with N, S z , J z ), with rrii — 0,1,2,3, 
mi > ?Ti2 > ... > rrif, and by filling the quanta m, of each partition in the boxes of a Young 
tableaux according to the following rules. 

i) The quanta must not increase when moving from left to right in a row or when moving down in 
a column. 

ii) The quanta referring to spin up and spin down states (m = 1, 2) must not appear more than 
once in a row. 

iii) The quanta referring to doubly occupied states or to empty states (m = 3, 0) must not appear 
more than once in a column. 

To pass from Young tableaux to states one usually apply standard Young symmetrizer and 
antisymmetrizer operators Q. Here, however, due to the fermion realization of the spin algebra, 
one must include minus signs in the permutations which involve the interchange of two fermions 
to account for the Pauli exclusion principle. Moreover, an additional (— l)- 7 sign is required for 
the permutations which involve doubly occupied states at site j to balance the alternating sign 
of the pseudospin algebra (without this sign the states would be eigenstates of S,S Z ,J Z , but not 
necessarily of J). To select among the filled Young tableaux the ones corresponding to highest 
weight vectors of 5*0(4), (actually highest weights for spin and lowest weights for pseudospin) we 
use the following criterion. Let us call the change of a spin down into a spin up a 1 — 2 flip and 
the change of a double occupied state into an empty state a 3 — shift. We note that S + just 
performs 1 — 2 flips while J~ performs 3 — shifts. It is clear that the filled tableaux which become 
unconsistent with the filling rules after 1-2 flips and 3-0 shifts arc annihilated by both S + and 
J~ i.e. they are highest-lowest weight vectors of the SO(4) algebra. On filled tableaux then, we 
perform 1-2 flips (i.e. turn a 1 into a 2 in all possible manners) to select the highest weights vectors 
of the spin SU(2), and 3-0 shifts to select the lowest weights of the pseudospin algebra. Note that 
according to our rules, double occupied states must be in the first row of a tableaux in consecutive 
order from the first box on the left. One can shift out a double occupied state of the first row 
without changing the spin symmetries, by moving the quanta one unit to the left and inserting a 
zero in the end box of the row, or by moving the quanta in the first column one unit upward and 
inserting a zero in the end box of the column. The tableaux that become unconsistent with the 
filling rules under both 1 — 2 flips and 3 — shiftings are the highest-lowest weight vectors of the 
SO (4) algebra. We remark however, that by 1 — 2 flips and 3 — shifts two different tableaux can be 
related to the same tableau of higher (lower) spin (pseudospin). In this case, linear combinations 
of the corresponding states must be taken. Having constructed the SO (4) invariant states we 
come to the problem of projecting them on a particular subgroup, say G, of Sf (we recall that any 
discrete group is a subgroup of Sf). To this end we denote by D(R),R £ Sf the irreps of Sf. A 
representation of G is readily obtained by selecting among the matrices D(R) those corresponding 
to elements of G. These representations however are in general reducible i.e. they can be expressed 
in terms of irreps D^> of G as D{R) — c u D^ v \R) with c„ non negative integers counting the 
number of times appears in D. By denoting with gi the number of elements in the class K% of 
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G and with g the order of this group, one easily express the integers c„ in terms of the characters 
X, X^"' °f respectively 5/ and G as 

Cv = -Y^9i~Ki ) *Xi- (5) 

^ i 

This gives the splitting of the irreps of 5/ into the irreps of G (here * denotes complex conjugation). 
The eigenfunctions i\) of S 2 , J 2 , corresponding to the above highest-lowest weight states, are then 
projected on the v— th irrep of G by using the projection operator pM 

= P("ty = ^ ^ x M *(iZ) f/R • ^ (6) 

where the sum is over all the elements R of the Sf subgroup (n v is the dimension of the z/-th irrep 
of G, x^"'{R) the corresponding characters and Ur the operator associated to the group element 
R) 0. By taking G to be the subgroup corresponding to the lattice translations T n we get the 
simultaneous eigenfunctions of S 2 , J 2 , T n with respect to which the Hamiltonian acquires a block 
diagonal form. To illustrate the method let us consider the case of the Hubbard model on a ID 
periodic chain with / = 4. This case corresponds to take the cyclic subgroup C4, of S4 which is 
an abelian group with ID irreps denoted by A, B, E\, E%. The irreps E\, E2 are one the complex 
conjugate of the other thus they are associated to the same energy level. We can take advantage 
of this accidental degeneracy by considering E\ , E2 equivalent to a single irrep E of dimension 
two (note that this degeneracy is connected with the time reversal invariance of the Hubbard 
hamiltonian). In Table 1 we have reported the filled tableaux associated to the highest- lowest 
weight vectors of the spin and pseudo-spin algebra constructed with our rules, for even values of 
the fillings (the odds fillings follow similarly and will be omitted for brevity). In this table we 
have also reported the splitting of the tableaux under the translation group. Note that, due to 
the presence of the (—1)-' sign in the Young operators, the even and odd representations of C4 are 
interchanged in the tableaux containing odd numbers of double occupied states. 

From Table 1 we see that the N=0 sector contains only one state (the vacuum) with eigenvalue 
Eq = 0. This state is a singlet under S and a quintet under J. The application of the pseudospin 
rising operator to it, generates an eigenstate of H of the N = 2 sector with S = 0, J = 2 but with 
J z = 1. The corresponding eigenvalue is E2 = U (the state does not appear in our table since it is 
not a lowest weight). By further applications of J + to the vacuum on gets all its descendents to 
higher (even) fillings 

N 

E N =E + -U, N — 2,4, 2f. (7) 

This is a general feature of the 50(4) diagonalization i.e. with the rising operator J + we can map 
all the spectrum computed at lower fillings N to the higher fillings TV + 2n, n — 1,2, .../ — N/2. 
From Table 1 we also see that the N — 2, S — 1 sector decomposes into two 1*1 blocks ( A, 
B states) and one 2*2 block (E states) for a total of 18 states, while the N — 2, S = sector 
decomposes into three blocks, two are 2*2 (B, E states) and one is 3 * 3 (A states) for a total of 9 
states. The dimension of the N = 2 sector is then 18 + 9 + 1 = 28 (note the addition of the state 
J = 2, J z = 1, coming from the N = sector). Similarly we see that the TV = 4 sector contains 
10 states with S = 0, 9 states with 5 = 1 and one state with 5 = 2, for a total of 42 states. If 
one adds these states to the ones coming from lower (even) fillings one gets the total dimension of 
the N = 4 sector as 42 + 18 + 9 + 1 = 70. These numbers just coincide with the ones obtained 
from the usual formula d,N = (// ) for N electrons on 2f sites. By rising these states to higher 
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fillings we get the total dimension of the Hilbert subspace corresponding to the even fillings as 
1*5 + 27*3 + 42 = 128. One checks that this number is just the sum of all the tableaux in Table 
1 taken with their multiplicities (i.e. the product of the tableaux dimension times the spin and 
pseudospin degeneracies) . This checks the completeness of our basis (the same analysis performed 
on odd fillings gives others 128 eigenvalues for a total of 4 4 states). In Tables 2a, 2b we have 
reported the SO (4) block decomposition of the hamiltonian matrix corresponding to the states 
listed in Table 1. Note that the eigenvalues at higher (even) fillings are obtained by diagonalizing 
the same blocks of Tables 2 but with an U added to the diagonal elements for each application 
of J + (the blocks at half filling obviously do not have descendents since J = 0). As expected, 
the ground state at half filling is a singlet with S = J = M. It is of interest to note that this 
state has momentum k = 7r (B state) and comes from the splitting of the same tableaux which 
characterizes the ground state of the 5/ invariant Hubbard model (see ||). One can easily 
compare these results with direct numerical diagonalizations in the full Hilbert space. Thus, for 
example, for N = 4 (di = 70) one finds that at t = .75, U = 1.5 the ground state is a singlet with 
energy E = —2.12132. This is just the same value obtained from the 5 = J = 0, B block of Table 
2a. 

Using the 5*0(4) symmetry we have also investigated a chain with / = 6 sites. Here for brevity 
we report only on the ground state at half filling (more details will be given elsewhere) . In this case 
we find that the ground state is obtained by diagonalizing a 14 x 14 band matrix corresponding 
to the 5 = J = 0, k = 7r block reported in Table 3. In this table the diagonal elements of the 
block are at the bottom while the elements above the diagonal (moving along columns), in the 
rows above. It is remarkable that the width of the band is just equal to the filling (this feature is 
true also for other fillings N < 6 Q). One can check that the eigenvalues of the 50(4) block in 
Table 3 are just the same as those obtained by diagonalizing the 400 * 400 matrix of the S z = 
space (for example, the ground state energy at t = 1, U = 1 in both cases is E = —6.60116). 

In closing this paper we shall estimate how the size of the 50(4) blocks will grow with /. This 
can be done by counting the number g?/,s,j of highest-lowest vectors of 50(4) for a fixed J, 5. We 
find that 

d f ,s,j - [(L) + (L- 2 )] - (L-x) [(U) + (L-i)] • w 

where c± = %± (5+ J) (the derivation of this formula will be given elsewhere |^]). One readily sees 
that Eq.(||) reproduces the correct number of states associated to the filled tableaux reported in 
Table 1. One also checks that these states, taken with their multiplicity, reproduce the full Hilbert 
space i.e. 

//2 f/2 

5(5 +1)J(J + 1)^ = 4/. (9) 

S=0 ,7=0 

By assuming an equal splitting of the 5/ representations into the irreps of the translation group, 
we estimate the dimension of the blocks as df,s,j/ f- Thus, for a chain with 10 sites the block 
characterizing the ground state at half filling is 560 x 560 while for a chain of 20 sites is 5*10 7 x5*10 7 . 
This shows that the 50(4) block diagonalization allows to reduce of about two order of magnitude 
the sizes of the matrices constructed by using just the S z and the translational symmetry (this 

last being j (/y T J (jv^))- Moreover, it is likely that (at least in the ID case) the band structure 
observed at / = 6 exists also at higher values of /. The possibility that the 50(4) diagonalization 
leads to band matrices for arbitrary values of / and N is presently under investigation. 
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Table Captions 



Table 1 

Decomposition of the filled Young tableaux corresponding to highest-lowest weight vectors of 50 (4) 
for / = 4 and N = 0, 2, 4. The splitting in terms of the irreps A, B, E, of the group C\ are shown. 

Table 2a 

50(4) decomposition of the Hamiltonian matrix corresponding to the highest weight vectors of 
the N — 2 states in Table 1. 

Table 2b 

50(4) decomposition of the Hamiltonian matrix corresponding to the highest weight vectors of 
the N — 4 states in Table 1. 

Table 3 

Band structure of the 50(4) block ( J = S = 0, k = ir) characterizing the ground state of the 
Hubbard chain with six sites. The diagonal elements of the block are reported in the row at the 
bottom while the elements above the diagonal (moving along columns) in the rows above. 
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Table 2a 
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Table 2b 
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TABLE 3 
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